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Abstract. We combine near-infrared photometry obtained with the VLT/ISAAC instrument and archival 
HST/WFPC2 optical images to determine VIK magnitudes and colours of globular clusters in two early- type 
galaxies, NGC 3115 and NGC 4365. The combination of near-IR and optical photometry provides a way to lift 
the age-metallicity degeneracy. For NGC 3115, the globular cluster colours reveal two major sub-populations, 
consistent with previous studies. By comparing the V — I, V — K colours of the NGC 3115 globular clusters with 
stellar populations models, we find that the colour difference between the two <; 10 Gyr old major sub-populations 
is primarily due to a difference in metallicity. We find A[Fe/H]= 1.0 ± 0.3 dex and the blue (metal-poor) and red 
(metal-rich) globular cluster sub-populations being coeval within 3 Gyr. In contrast to the NGC 3115 globular 
clusters, the globular cluster system in NGC 4365 exhibits a more complex age and metallicity structure. We 
find a significant population of intermediate- age very metal-rich globular clusters along with an old population 
of both metal-rich and metal-poor clusters. Specifically, we observe a large population of globular clusters with 
red V — K colours but intermediate V — I colours, for which all current stellar population models give ages and 
metallicities in the range ~ 2 — 8 Gyr and ~ 0.5Zq — 3Zq, respectively. After 10 Gyr of passive evolution, the 
intermediate-age globular clusters in NGC 4365 will have colours which are consistent with the very metal-rich 
population of globular clusters in giant elliptical galaxies, such as M87. Our results for both globular cluster sys- 
tems are consistent with previous age and metallicity studies of the diffuse galactic light. In addition to the major 
globular cluster populations in NGC 3115 and NGC 4365 we report on the detection of objects with extremely 
red colours (V — K <; 3.8 mag), whose nature could not ultimately be revealed with the present data. 
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Globular cluster systems are useful tracers of galaxy evo- 
lution. They consist of bright clusters, which are made of 
stars that share the same age and chemical composition. 
Globular clusters can form during major merger events in 



galaxies (e.g. Holtzman et al., 1992; Whitmore et al., 1993 



Whitmorc fc Schwcizcr, 1995 ; Schweizer et al., 1996| , etc!) 
but we also observe the formation of massive star clus- 
tcr in galaxies with a moderate star formation rate (e.g 
jO'Connell et al., 1994|; |Barth et al., 1995j jO'Conncll et al. 
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et al., 2000|, etc 
cluster 



1995; Brandl et al., 1996: Larsen & Richtler, 1999: Hunter 



In recent years the study of globular 



systems revealed the presence of globular cluster 



Table 1. Basic information on observed galaxies. The ref- 
erences are (1): |dc Vaucoulcurs et al. (199l[ ), (2) : |Schlcgc 



sub-po Dulations which must have formed in multiple for- 



Zepf. 



Harris 



mation epochs and/or mechanisms (see Ashman 
19981 ; |Kissler-Patig, 200C| ; |van den Bergh, 2000 
2001, far reviews). 

In practice, for distant galaxies, we observe integrated 
properties of their globular clusters, and would like to de- 
duce their physical properties, in particular their age and 
metallicity. Unfortunately, even in systems where redden- 
ing can be taken as uniform, the optically observed pa- 
rameters suffer from the well known age-metallicity de- 
gener acy. Spectroscopy can overcome this infamous prob- 
lem ([Jones fe Worthey, 1995| ; |Worthey fe Ottaviani, 1997| ; 



et al. (1998| ), ( 3): [Buta fc Williams (1995| ), (4): |Frogel et 
al. (1978| ). (5): |Tonry et al. (2001| ) 



Vazdckis fe Arimoto, 199£ ) to a large extent, but is, how- 
ever, very time consuming to perform for hundreds of ob- 
jects. Photometry still represents the most efficient way 
to study an entire globular cluster system. 

There have been several attempts to solve this degen- 
eracy for globular cluster systems in early-type galaxies 
which employed optical colours only (e.g. Kissler-Patig et 



Parameter 


NGC 3115 


NGC 4365 


Ref. 


type 


SO 


E3 


(1) 


RA (J2000) 


lOh 05m 14s 


12h 24m 28s 


(1) 


DEC (J2000) 


-07° 43' 07" 


+07° 19' 03" 


(1) 


I 


247.78° 


283.80° 


(1) 


b 


36.78° 


69.18° 


(1) 


Bt.o 


9.74 


10.49 


(1) 


Eb-v 


0.047 


0.021 


(2) 


(B - V)o 


0.94 


0.95 


(1) 


(V - I) eS , 


1.25 


1.25 


(3) 


(V - K) eB , 


3.30 ±0.02 


3.29 ±0.1 


(4) 


(m - M) v 


29.93 ±0.09 


31.55 ±0.17 


(5) 


M v 


-21.13 


-22.01 


(1),(5) 



al., 199^; |Whitmnrr rt al , 1 997[, [Kisslrr-Patig et al , 1008 



Kundu et al., 199E; Puzia et al., 1999) with various de 



the sensitivity to metallicity differences by a fa ctor of two 
or more (see Paper I, and SSP models of e.g. Bruzual fc 



Chariot 2000|) 



grees of success. Previous studies used the optical colour 
difference and the turn- over magnitude difference between 
the major globular cluster sub-populations to derive an 
age and metallicity difference. This method suffers from 
two main problems: 1) It is based on the assumption that 

the mqss fimptinns arp similar fnr hntn siih-pnpnlatinns 



In this paper we study the globular cluster systems 
of NGC 3115 and NGC 4365. NGC 3115 is an isolated 
galaxy, located at the very tip of the southern extension 
of the Leo Group, with just one si gnificant accompany ing 
nucleated dwarf elliptical galaxy ( Puzia et al., 2000a ) . It 
features a bimodal optical colour distribution of globular 



clusters (Kundu fc Whitmorc, 1998 ; Gebhardt fc Kissler- 



so thatuxumnosity differences directly rpfWt. Hiffwmcpg in 
the mass-to-light ratios of the individual populations. 2) 
The peak colour and turn-over magnitude differences arc 
small and somewhat model dependent. 

Alternatively one could derive mean age and metal- 
licity differences between globular cluster sub-populations 

from rplnu.r-rnlmir diagrams wnirh Hops not Hpppnn 1 on 



Patig, 1999| ; [Larsen et al., 200l|) 



NGC 4365 is a cluster giant elliptical galaxy located 
at the outer edge of the Virgo Cluster. Optical photom- 
etry of the globular cluster system revealed a broad but 
single-pea k colour distr i bution with little evidence for bi- 
modality ([Forbes 1996|; |Gebhardt fc Kissler-Patig, 1999 



the glnlhnla.r-p.lnRt.fr mass fimrtions 



Larsen et al., 200l| ; [Kundu fc Whitmorc, 2001a|) . purma fc 



Bender ( 1995 , see also Davies et al. 2001 ) detected a dc 



The combination of optical and near-infrared (near-IR) 
photometry can largely reduce the age-metallicity degen- 
eracy. This is because optical to near-IR indices like V — K 
are very sensitive to metallicity, but only have modest age 



sensiti vity ( Kissler-Patig, 200C ; Puzia et al., 2001 ; Kisslcr 



Patig £ t al., 2002, hereafter Paper I). Physically, this tech- 
nique works because the U-band samples mainly the light 
of stars near the turn-off, while the i^-band is most sen- 
sitive to cooler stars on the giant-branch in old stellar 
populations (see e.g. Yi et al., 2001). While the turn-off 



is mostly affected by age, the giant branch is primarily 
sensitive to metallicity (e.g. Saviane et al., 2000| ). Thus, 
V — I and V — K are affected similarly by age, but V — K 
is much more sensitive to metallicity (see Table [I] below) . 
This allows the age-metallicity degeneracy to be lifted by 
determining the location of a simple stellar population in 
a plot of V - I vs. V - K. 

Such diagrams are powerful tools to study the age and 
the metallicity of globular-cluster populations. Including 
near-IR passbands to our optical data we can increase 



coupled, counter-rotating core, which consists of a younger 
and more metal-rich stellar population than the rest of the 
galaxy. Such decoupled cores are believed to be created in 
major merger events. All relevant galaxy properties are 
summarized in Table [l]. 

The major goal of this paper is to derive ages and 
metallicities of globular clusters from a comparison of 
optical/near-IR colours with several SSP models. Putting 
this results into a larger context, including previous find- 
ings from galactic integrated light studies, we can further 
constrain formation scenarios for both galaxies. 

The outline of this paper reads as follows. In Section 
2 we describe the data reduction and calibration of 
our new near-IR VLT/ISAAC and the archival optical 
HST/WFPC2 data. The main results are presented in 
Section 3 followed by a discussion in Section 4. All ma- 
jor findings are summarized in Section 5. 

2. The Photometric Data 
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2.1. New VLT /ISAAC Near-Infrared Data 
2.1.1. Basic Reduction 

Deep images in were obtained in service mode (pro- 
gram 63.N-0287) with the Near-Infrared Spectrometer 
And Array Camera (ISAAC, see [Moorwood et al., 1998| ) 
attached to the Unit Telescope 1 (Antu) of the European 
Southern Observatory's Very Large Telescope (VLT). 
ISAAC is equipped with a Rockwell Hawaii 1024x1024 
pixels Hg:Cd:Te array. The pixel scale is 0.147"/pixel and 
the field-of-view 2.5' x 2.5' on the sky. In the following we 
will always refer to the K s filter as K. The ISAAC field 
was oriented such as to match the WFPC2 field of view. 

The NGC 3115 data were obtained during the nights 
of April 4th to 10th 1999 while the NGC 4365 data were 
taken during the nights of April 9th and 10th and ad- 
ditionally on the night of June 2nd 1999. The exposures 
were split into O — SSSSS — O sequences, where S was a 
dithered sky exposure comprising of 5 single exposures and 
O was a single on-target observation. Each sky image was 
exposed with a detector integration time (DIT) of 10 sec- 
onds and two such images were averaged in the read-out 
electronics (NDIT=2). Thus, the total integration time 
for a sky image is 20 seconds. An object image has a 
total integration time (DIT x NDIT) of 10 x 10 = 100 
seconds. Since the night-sky luminosity changes on short 
time scales the crucial part in each near-IR data reduction 
is the proper treatment of sky subtraction. Our observa- 
tion pattern allowed us to create mean sky images, which 
have been cleaned of stars prior to combination to a mas- 
tersky image. The data of each night have been processed 
individually. From each raw object image the mastersky 
was subtracted giving a sky/bias-subtracted object im- 
age. Those sky/bias-subtracted images were divided by a 
normalized masterskyflat, which is provided by the ESO 
quality control group. Eventually, all images were aligned 
using 15-20 objects common to all frames and averaged 
to give the final cleaned object image. The FWHM of 
the stellar PSF in the final K band image is ~0.6 " for 
NGC 3115 and NGC 4365. The total exposure time for 
NGC 3115 is 15500 sec, and 9500 sec for NGC 4365. 



2.1.2. Photometric Calibration 

For the photometric calibration of the NGC 3115 and 
NGC 4365 data set 11 and 9 near-IR standard stars have 
been taken throughout the nights ( Pcrsson et al., 1998| ). 
Each standard star was imaged five times, each exposure 
having the star centered at a different position on the chip. 
Standard-star photometry was performed with the pho- 



1 The K s (K short) filter is a slightly modified A'-band filter 
with the wavelength range 2.03 — 2.30/^m. Contrary to the K 
filter the transmission curve drops before the CO-absorption 
(at ~ 2.4/^m) band making the K s filter insensitive to CO 
abundance variations. Moreover, the sky background is reduced 
relative to the standard if-band due to the short-wavelength 
cut-off. 



tometry tool SExtractor v2.1.6 (Bcrtin & Arnouts, 1996) 
in an aperture of 30 pixels of diameter (7.4 x seeing). The 
corrections to an infinite-diameter aperture were found 
to be negligible in a subsequent curve-of-growth analysis. 
After applying a K — a clipping to all the measurements 
outliers have been deleted leaving 40 data points for the 
calibration of each data set. Due to the lack of sufficient 
spread in airmass among the standard-star measurements, 
we use the airmass term determined in the standard cali- 
bration of the ISAAC instrument for April 1999. Assuming 
that the colour terms are negligible in the near-infraredQ 
we find the following calibration relations for the photo- 
metric nights serving as reference 

K 3n5 = fc inst + 23.56(±0.018) - 0.05(±0.005)x (1) 
K 4 365 = fcinst + 23.75(±0.015)-0.05(±0.005)x (2) 



where -ftT ga iaxy is the calibrated magnitude, fci ns t is the in- 
strumental magnitude normalised to 1 second, and x the 
airmass. The error of the zero points includes photometric 
errors of each single standard star measurement and the 
errors of the curve-of-growth analysis. The error of the 
airmass term is a estimate from the variations in airmass 
of all single exposures. 

After correcting all nights to a common zero point (ob- 
tained from a reference photometric night) we traced the 
calibrated magnitude of two isolated bright stars over all 
nights. We find an upper- limit 1-er scatter being a few 
hundredths of a magnitude. Thus, we conservatively esti- 
mate the true photometric uncertainty, mainly driven by 
the strongly varying sky background, to be of the order of 
< 0.03 mag. 

Finally, all magnitudes were corrected for Galactic 
foreground reddening using the reddening values of Table 
|j] and the extinction curves of pardelli et al. (1989| ). The 
corrections for NGC 3115 and NGC 4365 are A K = 0.017 
mag and Ak = 0.008 mag. 



2.2. Archival HST/WFPC2 Optical Data 

NGC 3115 was observed with HST + WFPC2 under pro- 
gram GO. 5512, with the PC centered on the nucleus of 
the galaxy. The total exposure times of the combined im- 
ages are 3 x 350 sec in F555W, and 3 x 350 sec in F814W. 
NGC 4365 has been imaged with WFPC2 under program 
GO. 5920 in F555W and F814W filters with 2200 sec and 
2300 sec of total exposure time, respectively. The HST 
images were reduced and calibrated following for most 
parts the procedure as described in Puzia et al. (1999). 



Briefly recapitulating, all sub-exposures were matched by 



2 The theoretical value for the J — K s colour term was deter- 
mined by the data quality control group of the Paranal obser- 
vatory being close to zero (see ISAAC Data Reduction Guide 
1.5). Since no J-band data of NGC 3115 and NGC 4365 are 
available and the optical and near-IR photometric systems of 
Johnson and Persson do not overlap, i.e. they have no stars in 
common, it is not possible to include a V — K colour term into 
the calibration. 
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Table 2. Aperture corrections for WFPC2 photometry. 
The given terms refer to the correction from an 8/4- 
pixel-diameter aperture (PC/WF chips) to the standard 
Holtzman 0.5" aperture. # indicates the number of ob- 
jects which were used to determine the correction. 







PCI 


WF2 


WF3 


WF4 


<WF) 


NGC 3115 


V 


-0.152 


-0.296 


-0.385 


-0.323 


-0.323 




I 


-0.236 


-0.370 


-0.412 


-0.330 


-0.362 




# 


6 


25 


11 


24 




NGC 4365 


V 


-0.065 


-0.134 


-0.159 


-0.127 


-0.141 




I 


-0.133 


-0.205 


-0.242 


-0.219 


-0.225 




# 


28 


66 


39 


67 





whole-pixel shifts and stacked using a cosmic-ray cleaning 
routine (crrej) within IRAFQ. The object finding was per- 
formed on a master image which was created from all avail- 
able exposures of galaxy to allow reliable cosmic-ray clean- 
ing. None of the images suffers from crowding. Therefore, 
all magnitudes were measured with the SExtractor tool us- 
ing a 8 and 4-pixel-diameter aperture for the WF and PC 
chip, respectively. Aperture corrections to the Hol tzman 
standard aperture (0.5", see Holtzman et al., 1995| ) were 
determined in a curve-of-growth analysis. Only objects 
within the expected colour range for globular clusters were 
used (see Sect. 2.3). The correction terms can be found in 
Table |[ Instrumental magnitudes were then transformed 
to the Johnson V and / magnitudes according to the pre- 
scription given in Holtzman et al. 



All magnitudes were 
reddening corrected using the reddening values as listed 
in Tab le |l|. T ogether with the extinction curves of Cardclli 
et al. ( |L989| ) we obtain A v = 0.157 and Aj = 0.092 in the 
direction of NGC 3115 and A v = 0.070 and Aj = 0.041 
for NGC 4365. 

In the case of NGC 3115 excellent agreement wa s found 
between our data set and the previous studies of Kundu 
& Whi tmore (1998| ) and |Gebhardt fc Kissler-Patig (1999|) 
The systematic offset between the Kundu & Whitmore 
photometry and ours is 0.04 ± 0.11 mag in Vpsssw an( i 
—0.01 ± 0.11 in If814W- We find also good agreement 



for NGC 4365 between the photometry of Gebhardt & 
Kissler-Patig (1999) and our data (but see also Puzia et 
al., 2000b|). In Vfsssw we measure an offset 0.04 ± 0.07 
mag and 0.07 ± 0.05 mag for the /fsi4W filter. 



2.3. Selection of Globular Cluster Candidates 

For object selection the near-infrared ISAAC data were 
combined with the archival optical HST data. Note that 
the WFPC2 field size of 2.66' x 2.66' is well matched to 
the ISAAC field of view (2.5'x2.5'). After transforma- 
tion of HST coordinates to the ISAAC frame we re- 



3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Fig. 1. Selection parameter as a function of Vfsssw and 
Ifsuw magnitude. Plotted here are the full width at half 
maximum (FWHM) of the point-spread function (PSF) 
and the photometric error, as derived on HST images. 
Data points which were rejected during the selection are 
marked as open circles. Dotted lines indicated the selec- 
tion boundaries. 



quired the matches being better than 0.74" (~ seeing). 
Before selecting for likely GCs, the matched lists con- 
tain 110 and 186 objects with VIK photometry for 
NGC 3115 and NGC 4365, respectively. To get rid of ex- 
tended sources and noisy detections, we selected GC with 
cuts in FWHM\vfpc2 and photometric error, as measured 
on the optical HST data. Figure [l] illustrates the con- 
straints which define our globular cluster candidates. For 
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NGC 3115, all objects with a FWHM WFPC 2 larger than 
0.4" and for NGC 4365 larger than 0.2" were rejected. 
Both FWHM-WFPC2 values roughly correspond to 20 pc at 
the distance of 9.7 Mpc and 20.4 Mpc for NGC 3115 and 



NGC 4365, respectively (see fTonry et al, 200l| ). To deter- 
mine the real sizes of clusters which are rejected by the 
FWHM\vfpC2 cut we convolve the stellar PSF in each im- 



age with a series of King profiles (King, 1962) with a mean 
concentration parameter c = \og{r t /r c ) = 1.55, which was 
determined fr om the Milky Way globular cluster system 
( Harris, 1996 ) . The tests yield a maximum allowed core ra- 
dius r c w 5.6 pc and 4.9 pc, for NGC 3115 and NGC 4365 
globular clusters, to be retained in the sample. An in- 
spection of the core radii of Milky Way globular clusters 
shows that 82% and 79% have core radii smaller than the 
two respective maximum allowed sizes mentioned above. 
Unless the size distribution of NGC 3115 and NGC 4365 
globular clusters differs significantly, which seems unlikely 
given the similarity in sizes observed for many other early- 
type galaxies (e.g. Kundu fc Whitmorc, 199"§| ; Larscn et 
al., 20(0), we will miss only the extremely large globular 
clusters by applying the FWHMwfpc2 cuts. 

Furthermore, only objects with photometric errors in 
the HST photometry of less than 0.1 mag in V and I were 
retained. After the selection procedure the NGC 3115 list 
contains 87 candidates while the NGC 4365 list comprises 
136 globular cluster candidates^]. 

Background contamination is unlikely to make a sig- 



nificant contribution to the selected samples. Puzia et al 
(1999|)]have shown that with similar selection criteria the 
background contamination is less than 5%. 

2.4. Completeness 

The HST data for NGC 3115 is complete down to 
V F555W £ 26 and I F &uw £ 25. For NGC 4365 the 
photometry reaches Vf555W & 26.5 and Ifsuw 5s 25.5. 
Clearly, the limiting photometric passband is the if filter 
for both our galaxies and it will define the completeness of 
the data. Artificial star experiments with 1000 cycles, each 
including 100 artificial objects into the if images, yield a 
limiting magnitude (50% completeness) of if 5s 20.0 and 
if < 20.25 for NGC 3115 and NGC 4365, respectively. 

Despite the much longer exposure time for NGC 3115 
(15500 sec) the final if -band limiting magnitude is about 
the same as for NGC 4365 (9500 sec). Purely from the 
ratio of exposure times one would expect the NGC 3115 
photometry being 0.53 mag deeper than for NGC 4365. 
However, we find the NGC 4365 photometry is 0.25 mag 
deeper which adds up to a total difference of ~ 0.78 mag 
or a factor ~ 1.4 in S/N. 

2.4.1. Surface Brightness Fluctuations 

This discrepancy cannot be explained by a higher noise 
due to a difference in surface brightness between the halos 



of NGC 3115 and NGC 4365, since they are measured to 
be the same within 0.1 mag/arcsec 2 . However, a factor 
contributing to the lower detection limit in the NGC 3115 
field arc the larger if -band surface brightness fluctuations 
(SBFs) ( [Tonry fc Schneider, 1988| ): 

In addition to the Poisson noise of the flux, each res- 
olutions element contains SBFs which are the Poissonian 
fluctuations of the sampled number of stars. If the reso- 
lution element samples n stars of the mean flux F then 
the fluctuation of the flux will be a F = Fy/n. Hence, 
putting a galaxy twice as far away from the observer will 
keep its surface brightness constant (one resolution ele- 
ment samples now 4 times more stars which have 4 times 
lower flux) but lower the SBFs by a factor 2. The SBF am- 
plitude is proportional to the distance o> oc D . Both 
NGC 3115 and NGC 4365 have roughly the same if -band 
surface brightness, /X3115 = 12.5 and ^4365 = 12.6 mag 
arcsec~ 2 , at radii including the majority of our globular 
cluster candidates. However, since NGC 3115 is roughly 
two times closer than NGC 4365 flTonry et al., 2001| ), the 
ratio of SBF amplitudes, op oc S/N, will be ~ 2 (compared 
to ~ 1.4 measured). Thus, relaxing the implicit assump- 
tion of equal stellar populations, the difference in the SBF 
amplitude can explain the difference in the limiting mag- 
nitudes. 



2.4.2. Colour Completeness 

Regardless of the limiting magnitude, our data samples 
are biased toward red globular clusters. That is, the red 
(metal-rich) sub-population is more complete than the 
blue (metal-poor) one at a given magnitude threshold (see 
Fig. ||). To quantify the biasing factor we assume that 
the globular cluster luminosity functions (GCLFs) of both 
sub-populations have the same shape (same mass func- 
tion) and the same dispersion (same slopes of the mass 
function). Both globular cluster sub-populations are as- 
sumed to have the same M/L. Previous studies of the glob- 
ular clusters systems of elliptical and SO galaxies suggest 
that this assumption is a reasonable initial approximation, 
although in detail there is some evidence for modest dif- 
ferences in the M/L between the metal-poor and metal- 
rich sub-populations of well-studied cluster systems (see 



Puzia et al., 1999; Kundu et al., 1999, for a discussion on 



The final lists are available in electronic form from THP. 



NGC 4472 and M 87). 

The turn-over magnitude of the NGC 3115 GCLF is 
Vto = 22.6 ± 0.2 for both blue and red cluster candidates 
( Larsen et al., 200l] ) . Concerning the mean V — K colour 
for the blue (2.3 ±0.1 mag, see Sect. [D] and Table g) and 
for the red clusters (3.0 ±0.1 mag) the expected turn-over 
magnitudes in if are if to, blue ~ 20.3 and if to, red ~ 19-6 
mag. According to the completeness limit of the NGC 3115 
data (ifii m w 20.0 mag) and assuming a GCLF dispersion 
cgclf = 1-3, we sample 41% of the blue and 62% of the 
red sub-population. The biasing factor for the NGC 3115 
globular cluster system is thus ~ 1.5. That is, we detect 
^1.5 times more red clusters than blue ones. In the case of 
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NGC 4365 we detect 13% of blue and 23% of red clusters 
assuming Vto = 24.37 ± 0.16 mag along with ctqclf = 
1.22 mag ( |Larsen et al, 200l| ) and V - K = 2.76 and 3.23 
mag for blue and red clusters, respectively. The biasing 
factor is 1.8. If we assume a single peak colour distribution 
with a mean V — K « 2.95 (see Sect. 3.1) and a turn-over 
magnitude independent of colour we sample 17% of the 
GCLF in NGC 4365. 

2.4.3. Spatial Completeness 

Both our globular cluster samples are spatially incomplete 
(see Fig. ||). Moreover, the different spatial distribution of 
blue and red globular clusters also influences the com- 
pleteness. In NGC 3115 we sample the inner 4 kpc of the 
globular cluster system. The spatial coverage of NGC 4365 
globular clusters is more complete. Our field of view covers 
clusters wit h galactocentric distances up to ~ 10 kpc. For 
NGC 3115, |Kundu fc Whitmorc (199S| ) report that the red 
globular clusters follow the galactic light profile while the 
blue clusters rather reside in the halo. Thus, independent 
of co mpleteness biases, both samples have to be consid- 



Table 3. KMM mean colours of the globular cluster sub- 
populations in NGC 3115 and NGC 4365. 







V 




K 


V 


_ ja 


NGC 3115 


blue 


2.31 


± 


0.06 


0.922 


±0.03 




red 


3.01 


± 


0.04 


1.153 


±0.03 


NGC 4365 


blue 


2.76 


± 


0.04 


0.981 


±0.03 




red 


3.23 


± 


0.04 


1.185 


±0.03 



values for the blue and red sub-populations adopted from 
Larsen et al. (200l|) . 



ered a a spatially biased in favour of red globular clusters. 
The spatial distribution of globular cluster candidates is 
shown in Figure §, along with the location of the WFPC2 
and ISAAC apertures. 

3. Results 



3.1. C olour Distributions and Colour Magnitude 
Diagrams 



We use the KMM code ( |Ashman et al., 1994| ) to test for 
bimodality of the colour distributions. The code calcu- 
lates the likelihood Pkmm whether or not a two-Gaussian 
fit is preferred over a single-Gaussian fit and determines 
the positions of the modes. Higher Pkmm values indicate 
that a single-peak distribution is more likely. In addition 
to the cuts in FWHMwfpC2 and photometric error (see 
Sect. |2~3| ), we select globular cluster candidates by their 
V — K colour in the range 1.7 < (V — K) < 3.7 and 
cut the samples at the limiting K magnitude (NGC 3115: 
K < 20.0, and NGC 4365: K < 20.25). According to 
simple stellar population models (see e.g. Fig. |) this 
colour range selects globular clusters with metallicities 
—2.5 ^[Fc/H]^ ±0.4 dex and with a wide range of 
ages. The input lists contain now 69 and 78 objects for 
NGC 3115 and NGC 4365, respectively. 

KMM returns for NGC 3115 a likelihood of P K mm = 
0.53 with peaks at V - K = 2.31 ± 0.06 and 3.01 ± 0.04. 
Former studies found a clear bimodality in the V—I colour 



distribution of globular clusters in NGC 3115. Kundu 



Whitmore (1998| ) find peaks at V — I = 0.96 and 1.17 and 
the most likely separati on of the two sub-populations be - 
ing at (V-J)dip = 1-06. |Gebhardt fc Kissler-Patig (1999|) , 
who first calculated the likelihood of a dip in a supposedly 



peaks^] at V - I = 0.92 ± 0.03 and 1.15 ± 0.03 and a dip 
likelihood P dip = 0.96. 

As for NGC 4365, KMM gives Pkmm = 0.94 with 
peaks at V-K = 2.76±0.04 and 3.23±0.04. In NGC 4365 
|Gcbhardt fe Kisslcr-Patig (1999| ) find a hint for bimodal- 
ity in their V — I colour distribution. The mean colour 
and the bimodality significance of their distribution are 
(V ~ I ) mean = 1.09 ± 0.02 and P dip = 0.76. [Larson ct al 
(2001P measure V — I = 0.98 ± 0.03 and 1.19 ± 0.03 as 



peak colours and a dip likelihood Pji p = 0.033. These re- 
sults show that there are some galaxies for which extant 
colour distributions do not give a clear result when statis- 
tical tests are applied to search for bimodality. No clear 
multi-modalit y is detectable in cither colour d i stribution 
(for V — I sec Gcbhardt fc Kissler-Patig, 1999 ; Larsen et 



al., 200 1| ; |Kundu fc Whitmore, 2001a|) . The globular clus- 



ters show a continuous colour distribution. 

The errors of V — K peak colours in our KMM analysis 
have been determined by varying the input parameter of 
the KMM code, which are the initial colour for both peaks, 
their covariance, and whether or not the two peaks have 
the same dispersion. All values are summarized in Table [s| 
Note, that the driving parameters for detecting bimodal- 
ities with KMM are the peak locations, the dispersions, 
and the sample size. With smaller samples the Pkmm pa- 
rameter becomes less meaningful ( Ashman et al., 1994 ). 
The bimodality test shows that only the colour distribu- 
tion of NGC 3115 is likely to be not unimodal. Contrary, 
NGC 4365 has a globular cluster colour distribution, which 
is well represented by a single Gaussian. See Figure || for 
(V — K) histograms. 

A metallicity calibration of V — K using Milky Way 
and M31 globular clusters with Eb~v < 0.27 is derived 
in Paper I 

[Fe/H] = -5.52(±0.26) ± 1.82(±0.11) -{V-K) 

with a rms of 0.29 dex. With this calibration we calculate 
metallicities from peak colours which were obtained in the 
KMM analysis. The globular cluster sub-populations in 
NGC 3115 have then metallicities [Fe/H]= -1.32 ±0.11 
and -0.04 ± 0.07 dex. Transforming the KMM peaks of 



-I) 



1.03±0.03 



bimodal colour distribution, give (V- 
with a dip likelihood Pdi p = 0.71. Larsen et al. (2001 ) find 



5 The errors of the colour peaks were kindly provided by 
S0ren Larsen. 
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Fig. 2. Field of view for the NGC 3115 (left panel) and NGC 4365 (right panel) data on the sky. The solid square shows 
the effective aperture of the ISAAC data. The HST aperture is indicated as a dashed polygon. Solid ellipses indicate 
the area encompassing half of the integrated light of the galaxy. For NGC 3115, we split the sample at V — I — 1.05 
and plot the blue sub-population as open triangles and the red clusters as filled squares. Stars mark the position of 
extremely-red objects in the colour range 3.7 < V — K < 4.3, which meet our completeness but not the size selection 
criteria (all extremely-red objects lie above the dashed line in Fig. ^). For NGC 4365, solid circles show the positions 
of globular clusters with intermediate ages, whereas the remainder globular cluster candidates are marked by open 
circles. 



the NGC 4365 system yields [Fe/H]= -0.50 ± 0.07 and 
0.36 ±0.07 dex. The median V — K colours of the globular 
cluster systems of NGC 3115 and NGC 4365 are equivalent 
to [Fe/H]= -0.15±0.07 and -0.33±0.05 dex, respectively. 
The errors include the uncertainty in the peak colours 
only. Note, that the calibration was established with Milky 
Way and M31 globular clusters, which have very few clus- 
ters at higher mctallicities (i.e. at [Fe/H]<; —0.4 dex). 
Consequently, all derived values with metallicities around 
solar and higher may suffer from uncertainties due to ex- 
trapolations. 

In Figure || (bottom panels) we show V vs. V — K 
colour-magnitude diagrams (CMDs). Together with the 
V — K histograms in Fig. || (top panels) the CMDs of 
NGC 3115 and NGC 4365 show no clear multi- modality 
due to the limited sample size. The dashed lines show the 
50% completeness limits and make the selective complete- 
ness obvious. 

The CMDs also show the data which were rejected by 
our object selection criteria (FWHMwfpc2 , photometric 
error, and V — K), which are shown as open circles. For 
NGC 3115, this includes two bright and blue likely fore- 
ground stars. Using the galactic stellar population model^] 
of Robin et al. (1996| ), we expect zero stars. In the case 



of NGC 4365 we find no blue contaminating objects, in 
agreement with the model predictions. 



3.2. Extremely Red Objects 

At the red end of the NGC 3115 and NGC 4365 CMD 
(Fig. |3|) we find a dozen objects with V — K s» 4.0 mag, 
which were rejected by the FWHM\vfpc2, completeness, 
and photometric-error cuts. These objects show no signifi- 
cant clustering and are evenly distributed on the sky. Some 
of them show an extended PSF and their FWHMwfpC2 
values exceed our FWHM-threshold making them likely 
background galaxies. The nature of the remaining very red 
point sources (indicated by squares in Fig. ^), not rejected 
by these selection criteria but only by the V — K colour 
cut, is less clear. We discuss in turn some possibilities of 
their nature: 

Foreground stars: In the colour range 3.5 Si V — K ^ 
4.5 mag and the magnitude range 22.0 ^ V & 24.5 mag we 
expect two stars in the ISAAC field of view (6.25 arcmin 2 ), 
according to the model of [Robin et al. (1996| ). The effec- 



6 The model predictions were calculated using the code at 
www.obs-besancon.fr/www/modele/modele_ang.html 



tive field of view is even smaller due to the incomplete 
overlap between the ISAAC and WFPC2 field. That is, 
significant stellar contamination at this red colours and 
faint magnitudes can be considered as unlikely. 

Fuzzy red globular clusters: In optical colours the 
extremely- red objects are clustered around V — I ~ 1.4 
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Fig. 3. V vs. V — K colour-magnitude diagrams for NGC 3115 (upper panel) and NGC 4365 (lower panel). The bottom 
of each panel shows the CMD of the associated globular cluster system. Filled symbols are selected GC while open 
symbols indicate data, which were rejected by our selection criteria (i.e., FWHMwfpC2, photometric error, and colour). 
The dashed line shows the 50% completeness limit of the photometry as constrained by the near-IR photometry. The 
errors include photometric errors only. The top of each panel shows the V — K histogram of all detected sources (solid 



histogram) and the selected GC sample (shaded histogram). A probability-density estimate (Silverman, 1986; Puzia 



et al., iiOOOa) of the selected unbinncd data is indicated (thick line) along with its la uncertainty limits (dotted lines). 
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Fig. 4. Relative object sizes as a function of V — K colour 
for NGC 3115 (top panel) and NGC 4365 (lower panel). 
Solid points show the data after application of the selec- 
tion criteria (FWHM W fpc2, photometric error, and U— K 
colour) while open circles show all data. Open squares 
show data with red colours which would pass the selection 
criteria if the V — K selection were dropped (i.e. point- 
like very red objects). Note, that there is no trend of size 
vs. colour see n in either galaxy, contrary to the results of 
Larsen et al. who find that red globular clusters in near- 
by elliptical galaxies are on average 20% smaller than blue 
ones. However, this correlation might be camouflaged by 
the colour bias of our sample. 



and 2.2 mag. Counterparts of the redder objects could 
be the recently detected population of extremely-red ob- 
jects around NGC 1023, with typical colours V — I w 2.2, 
which appear to be more extended than the other glob- 
ular cluster candidates ( |Larsen fc Brodic, 2000 ). Using 
the resolution power of HST, we measure the A(0.5 — 3) 
parameter (which is simply the magnitude difference be- 
tween the flux within two different aperture sizes; see e.g. 
Puzia et al., 1999) to study relative sizes of our globular 
cluster candidates. Figure [| shows the relative size distri- 
bution of globular cluster candidates and all other objects 
of our samples. The wide-field chips of WFPC2 on-board 
HST resolve objects down to 4.5 and 9.6 pc radius at the 
distance of NGC 3115 and NGC 4365, respectively. Red 
objects, indicated by an open square (in Fig. |j), were re- 
jected by the V — K colour selection only. 

Intermediate- age globular clusters: We point out that 
the presence of carbon-star dominated intermediate-age 
globul ar clusters, like they a re found in the LMC and 
SMC ( Olszewski et al., 1996 ), is unlikely to explain the 
extremely red colours. The strong IR flux of such AGB 
stars, yields red V — K colours but leaves the V — I colour 



relatively unchanged. The V — I colours of the extremely- 
red objects are too red to be matched by any simple stellar 
population at low redshift, regardless of its age and AGB 
contribution. 

Background galaxies: Unresolved background galaxies 
can, however, also mimic a stellar PSF. In the case of 
our very red objects, we might deal with the blue end of 
the known population of Extremely Red Objects (EROs), 
identified as red-shifted and/or obscured galaxies at z > 
0.85 (e.g. [Liu et al., 2000] paddi et al., 200l| ). The colour 
definition of EROs varies in the literature and is typically 
around R — K <; 4.5 — 6 mag. The surface density of 
EROs is a function of colour and limiting magnitude, and 
increases from 0.07 arcmin -2 for R — K > 6 mag to 0.49 
arcmin -2 for R—K > 5 mag for a sample with Ky lm < 19.0 



mag (Daddi et al., 2001). In our field of view, we expect 
0.4 and 3.1 EROs, respectively, down to K = 19 mag. Our 
objects have similar colours as the bluest EROs of which 
the distribution on the sky is inhomogeneous. Hence, the 
number of very red objects in our field of view (two at K ~ 
19 mag) is consistent with the surface density of EROs. 
The red end of the NGC 4365 CMD shows also hints of 
an extremely-red population of objects fainter than the 
photometric limit. All, except two sources in this colour 
regime are rejected by the FWHMwfpC2 selection and the 
completeness limit. 

In conclusion, the nature of our extremely-red objects 
cannot be definitely clarified from the photometric data 
only. Spectroscopy is necessary to distinguish between cool 
M-type foreground stars, obscured background galaxies, 
or (unobscured) extremely metal-rich young globular clus- 
ters. 



3.3. Colour-Colour Diagrams 

With the VIK data in hand optical/near- IR colour-colour 
diagrams, such as V — I vs. V — K, can be constructed. 

In the following we compare globular cluster colours 
to different simple stellar population (SSP) models and 
determine the mean age and metallicity of the major sub- 
jopulations. We use the models of Vazdckis (1999 ), Kurth 
2t al (1999| ), [Bruzual fc Chariot (2000| ), and |Maraston 
(200l| ). Absolute calibrations of SSP models are nowa- 
days still uncertain. However, the model-to-model differ- 
ences of relative values are small. This is particularly true 
for the age and metallicity ranges covered by Milky Way 
and M31 globular clusters which serve as model calibra- 
tors. Figure |5| demonstrates model-to-model variations for 
a stellar population of 12 Gyr with a metallicity range 
-1.7 <;[Fc/H]<; +0.4 dex. Also plo tted are the globular 



clusters of M31 and the Milky Way ([Harris, 1996^ [Barmby 
|ct al., 2000| ; |Aa7 onson fc Malkan, in prcp.p . The plot shows 
that at a fixed age, metallicity increases along the track to- 
wards redder colours. Note, that there is very few data for 
M31 or Galactic globular clusters in the solar-metallicity 
regime. Tracks of younger ages fall generally below (i.e., 
have bluer V — I colours than) older tracks and the colour 
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Fig. 5j Variationa of model prcdictiona for a 12 Gyr old 
simple btcllar population in the mctallioity range [Fc/H] = 



alone. The second result is that the two globular cluster 
sub-populations have both old ages with a small or a negli- 
gible age difference. Overall, the whole system seems to be 
dominated by old globular clusters. The colour difference 
is, therefore, largely due to a metallicity difference. 

To quantify the age and metallicity difference between 
both globular cluster sub-populations, we compare the dif- 
ference in their V — I and V — K colours to several SSP 
model predictions. We recall again that our data sample 
is biased towards red globular clusters, i.e., we detect 1.5 
times more red clusters in the NGC 3115 system than blue 
ones (see Sect. 2.4). The blue peak of the V — K colour 
distribution might, therefore, vary as one goes deeper with 
the K band photometry, adding more blue globular clus- 
ters to the sample. 

Figure (?] shows the mean colours of the two sub- 
populati ons ovcrplotted on SSP model grids. From the 
model of [Maraston (2001 ) we derive an age difference be- 
tween the two sub-populations of At = (tbiue - tred) = 2+4 
Gyr. The other models yi eld age differen ces of — 5 ± 6 
([Bruzual fc Chariot, 2000|) , -4+5 Gy r flKurth et aT 



-2.25 to +0.35 dex. The lines indicate the predictions 



of Maraston (2001) (solid), Kurth et al. (1999) (short 



dashed), |Bruzual fc Chariot (2000[ ) (dotted), and jvazdekis 
(1999|)](long-das hcd) . The metallicity on the iso-age line of 
Maraston (2001) is labeled and indicated by filled squares. 



1999J ), and -4 ± 5 Gyr ( |Vazdekis, 1999[ ). The mean age 
difference is (At) = — 2± 3 Gyr. The two globular cluster 
sub-populations are coeval within la. 

With the assumption of coeval sub-populations 
we calculate the metallicity difference which is the 
driving parameter for the V — K colour offset. 



Open triangles are globular clusters of M31 flBarmby et al 
2000|)~|vhi le open circles mark the colo urs of Milky Way 
globulars ( |Aaronson fc Malkan, in prep. ). The (V—I) data 
was taken from Harris (1996 ). 



Applying the model of Maraston (2001) yields A[Fe/H] = 
|[F e/Hlbiuc-[Fe/Hl rcd | = 1.48 ± 0.57 dex. The models 



of Bruzual fc Chariot, Kurth et al., and Vazdekis give 



differences between tracks increases with younger ages. 
For old ages the separation of tracks shrinks to values 
comparable with photometric errors which diminishes the 
ability to distinguish between ages. 

3.3.1. NGC 3115 

The colour-colour plot for NGC 3115 globular clusters in 
Figure [j| show selected data with their photometric er- 
rors. Large squares show the mean colours of the two sub- 
populations as determined in the KMM analysis. Their 
error bars indicate the error of the mean as defined by 
Cmean = a / yN . As a guidance for the eye, iso-age and 
iso-me tallicity lines taken from the Bruzual fc Chariot 
|(2000] ) SSP model are overplotted. Dotted lines show iso 
age tracks for 1 to 15 Gyr populations with metallicities 
-1.7 <[Fe/H]< +0.4 dex. 

Figure |^ reveals two main results. First, there is a 
broad spread in metallicity with evidence for bimodality in 
the NGC 3115 cluster system. The globular cluster popu- 
lation can be split at around V — I ~ 1.05 and V—K 2.6 
into two major sub-populations. In general, the colour- 
colour diagram helps to pick out the multi-modality which 
is not evident in V — I or V — K colour distributions 



1.04+0.30, 1.03+0.31, and 0.99+0.19, respectively. The 
weighted mean metallicity difference is A[Fe/H]= 1.0+0.3 
dex. The metallicity difference derived wit h the [Fe/H]- 
(V-K) relation from Paper I (see Sect. |+]) is 1.28 + 0.32 
dex and is consistent within the errors with the model 
predictions. 

3.3.2. NGC 4365 

The NGC 4365 colour-colour diagram is more complex 
than the one for the NGC 3115 globular cluster system. 

Dust: Can extinction by dust in NGC 4365 produce 
the colour sequence we observe? The reddening vector is 
plotted in the upper left corner of Figure ||. If extinction 
by dust in NGC 4365 would play a role at all, one would 
expect this to have more impact within the galaxy body 
than in the halo, which would cause the bulge population 
to be systematically redder. Correcting for any additional 
internal reddening would shift the clusters along the red- 
dening vector almost parallel to the iso-metallicity tracks 
towards younger ages (see Fig. [j]). However, there is no 
evidence for the existence of significant amounts of dust 
in NGC 4365 ( [Jura et al., 1987| ; poudfrooij et al., 1994| ). 

Intermediate- age clusters: The most striking feature 
in Figure ^ is a large population of globular clusters with 
colours that are red in V — K (~ 3.0) but intermediate 
in V - I (~ 1.0 - 1.2). According to SSP models this 
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Fig. 6. V — I vs. V — K colou r-colour diagrams for NGC 3115 and NGC 4365. All data are reddening corrected 
according to Sect j2. 1.2] and 2/2. For comparison, we show six isochrones for simple stellar populations of 15, 10, 5, 3, 
2, and 1 Gyr (dotted lines from red to blue V — I colours). The dashed lines represent iso-metallicity tracks changing 
from the lower left to the upper right from [Fe/H]= —1.7 to +0.4 dex. The SSP model used here was adopted from 
Bruzual fc Chariot (2000| ). For NGC 3115, the large squares mark the mean colour of the modes as determined by 
KMM code (see Sect. |5j). Arr ows in the upper left corner of each panel indicate the applied extinction correction and 
the direction of reddening. 
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Fig. 7. NGC 3115: V — J vs. — K colour-colour diagrams for major globular cluster sub-populations. The four 
panels show the mean colour of the two major globular cluster sub-populations in NGC 3115 with four different 
SSP model grid. In the upper left panel with the SSP model of Maraston (200l[ ) we labeled iso-age (solid lines) and 
iso-metallicity tracks (dashed lines) in the range 2 — 15 Gyr and [Fe/H]=— 2.25 to +0.35 dex. All other model grids 
have exactly the same iso-age lines. For better comparison we plotted the 12 Gyr iso-age line thicker. Triangles show 
the position of the labeled metallicity along the 12 Gyr age track. The SSP models of Bruzual & Chariot (2000); 
Kurth et ah (1999); Vazdekis (1999| ) have slightly different metallicity spacings. The are iso-metallicity tracks for 
[Fc/H]= +0.4, 0.0, —0.4, —0.7, —1.7, and —2.3 dex starting at red colours (upper right in each panel). 



clusters fall in the regime of ~ Zq — 2>Zq metallicities 
and intermediate ages of ~ 2 — 8 Gyr. In general, the 
red globular clusters appear on average younger than the 
blue sub-population.[] With a few exceptions all metal-rich 



7 We briefly review results of previous GCLF studies which 
were based on optical photometry only. From recent HST pho- 



([Fe/H]> -0.4 dex) clusters have V — K, V — I colours 
indicative of intermediate to young ages. We estimate this 
intermediate-age clusters to represent ~ 40 — 80% of our 
observed globular cluster sample. This limits are derived 
1) by assuming all clusters red of the [Fe/H]= —0.4 dex 



our findings. However, we note, that the V — I colour distri- 



tometry of 500 globular clusters 

difference in GCLF turn-over magnitudes ArriTO = 0.8 ± 0.4 
mag between the red and blue sub-population which they cut 
at V — I — 1.05 mag. With the red clusters being fainter, 
their result implies that the red clusters are on average older 
and more metal-rich than the blue sub-population, contrary to 



Larsen et al. (2001) find a bution (see Fig. 4 in Larsen et al. 2001 ) is very broad with 



no apparent bimodality. Since Larsen et al. 's V — I separation 



splits between blue and red clusters at roughly the colour of 
the intermediate-age globular cluster population, their GCLF 
analysis might be corrupted by the intermediate-age clusters 
that were not accounted for. 
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iso-metallicity line (see Fig. to be intermediate-age and 
2) by taking into account the existing spread in colours of 
globular clusters in Local Group galaxies and only count- 
ing objects with colours redder than all Milky Way and 
M31 globular clusters and with V — K <; 3.0 mag as 
intermediate-age clusters. 

SSP models predict that the V magnitude of a ~ 2 — 8 
Gyr globular cluster is ~ 1.7 to 0. 6 mag brighter than 



that of a 15 Gyr old cluster (e.g. [Bruzual fc Chariot 
2000|)~|Thus we should be able to detect bright globu- 
lar clusters in the CMD. The NGC 4365 CMD (lower 
panel in Fig. H) shows indeed a population of relatively 
bright globular clusters around 21 ^ V 22 mag and 
2.6 ^5 V — K^i 3.4 mag. Most of these objects, indeed, fall 
in the colour-colour diagram (see Fig. ^) into the regime 
of metal-rich intermediate-age populations with ages ^ 8 
Gyr and metallicities Z <; Z & . As a consistency check 
we, conversely, select clusters in the colour-colour diagram 
(Fig. with colours of metal-rich intermediate-age popu- 
lations (see colour definitions above) and look for their V 
magnitudes. They dominate the GCLF at brighter magni- 
tudes: At V ^ 22.0 mag ~ 80% of all objects in the sample 
are metal-rich and have intermediate ages. Unfortunately, 
the small size and the colour bias of our sample does not 
allow to establish statistically whether or not the GCLF 
of the metal-rich intermediate-age population is brighter 
than that of the remaining globular clusters. 

Comparison of ages with known young clusters: 
Clusters of pre-AGB ages (sS 100 Myr) would have not 
only different colours, but also considerably brighter mag- 
nitudes than the old population and would thus be easy to 
identify. A lower limit on the ages of the intermediate-age 
population can be set by comparing their colours with 
the data fo r very young (;$ 1 G yr) globular clusters in 
NGC 7252 flMaraston et al., 200l|) . All NGC 7252 clusters 
show too blue V — I colours compared with the bulk of 
intermediate-age clusters in NGC 4365. This implies that 
the latter are at least 1 Gyr old. 

A recent study of globular clusters in NGC 1316 re- 
vealed an intermediate-age (3 ±0.5 G yr), solar- metallicity 
globular cluster sub-population ( |Goudfrooij et al. 
2001a[jp|). The brightest and not-reddened intermediate- 



age cluster candidates have mean colours V — I ~ 1.05 
and V — K w 2.87. This colours, compared with the mean 
locus of globular clusters in NGC 4365 (cf. Fig. ||), indi- 
cate that some of the intermediate- age globular clusters 
in NGC 4365 have counterparts in NGC 1316 which are 
- 3 Gyr old. 

Spatial distribution: Interestingly enough, the young 
globular clusters with colours V — K 2.6 seem to form 
a 'bulge' or 'disk' population, as the distribution of those 
clusters is clearly elongated along the major axis of the 
galaxy body (cf. Fig. [^). The major axis of the disk 
aligns with the major axis of the galaxy body. In con- 
trast, blue globular clusters form a halo population around 
NGC 4365. Unfortunately, our field of view is concentrated 
to the south-east only and does not allow to check for con- 



sistency on the opposite side (north-west) of the major 



3.4. Aging the NGC 4365 Globular Cluster System 

In the following we let the entire globular cluster system 
in NGC 4365 age for 10 Gyr and study their evolution 
in the colour-colour diagram. We assume a passive colour 
evolution due to age only. Each single globular cluster is 
aged individually for 10 Gyr. 

Following Buzzoni (1995) a detailed inspection of each 
model shows that the parameter range 5 ^ t ^ 15 Gyr and 
—0.7 Ss[Fe/H]^S 0.4 dex (except for Maraston models, for 
which we use a smaller metallicity range —0.35 <^[Fe/H]<i 
0.35 dex) allows a simple parameterization of colours by 
age and metallicity. The parameterization can be written 
as a first approximation: 



AColour, = oiiAlog t + /3 4 A[Fe/H] 



(3) 



where the slopes a*, (3i are determined from the SSP mod- 
els and are given in Table ||. We use linear interpolation 
in each model grid and average the obtained ages and 
metallicities from all four models. For clusters with redder 
colours than the model parameter space (see Table ||) we 
use an extrapolation to high-metallicities, being aware of 
the insecure prediction power of SSP models in this range. 
Clusters which have bluer colours than the model param- 
eter space are dropped from the evolution due to apparent 
non-linearities in the grids beyond [Fe/H]^ —0.7 dex and 
an expected small colour evolution. Since we assume pas- 
sive evolution one can write 



A{V -I) = 0.230 log(l 
A(V-K) = 0.514 log(l 



lOGyr 
t 

lOGyr 



(4) 
(5) 



where t is the current age of each single globular cluster de- 
termined from the colour-colour diagram and A(V— I) and 
A(V — K) are the colour changes over a 10 Gyr lifetime. 
The mathematical form of these equations implies a rapid 
evolution in colour for intermediate-age globular cluster. 
As the clusters get older, the colour evolution slows down, 
i.e. old clusters are hardly affected in colour by the aging 
procedure. 

We emphasize again that both upper equations are 
only valid for stellar populations older than ~ 5 Gyr (see 
Table ||) and that, especially for younger ages, the rela- 
tions represent insecure extrapolations. Figure || shows the 
evolution of the colour-colour diagram and the according 
colour distributions. 

Only clusters which fall right to the dot-dashed line are 
considered in the aging process. This boundary is defined 
by the limitation of our parameter space (see Table [3]). 
Moreover, there are very few objects left of the line for 
which the models predict almost no colour evolution. The 
aging shows that the 10 Gyr older globular cluster system 
(solid dots) contains a significant fraction of clusters with 
Z > Zq metallicities. Clusters with extremely red colours 
V — K > 4, which correspond to metallicities Z <; IOZq, 
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Fig. 8. NGC 4365: Colour-colour diagram for the evolved globular cluster system. Solid circles show the evolved 
data while open circles indicate the colours of the observed data set. Since an average evolution calculated from four 
SSP models was adopted we plot 12 Gyr iso-age tracks for each model with mctallicitics —2.2 5 <[Fe/ H]< 0.35 dex 
The metallicity seque nce is indicated by solid s quares. We us ed models from |Maraston (2001 ) (solid), Kurth et al 



(1999|)1(short-dashed), [Bruzual fc Chariot (2000|) (dotted), and |Vazdekis (1999|) (long-dashed). All clusters right of the 
dot-dashed line, which is the limit of the applied linear interpolation, were aged. 



are present in the new colour-colour distribution. Note, 
however that these extreme colours are the predictions 
of the extrapolated parameter range and might change 
when detailed models of extremely metal-rich stellar pop- 
ulations are available. 

The aged V — I and V — K colour distributions (thick 
lines Fig. ||) show a very broad flat-top red peak. This 
broad peak suggests a wider metallicity dispersion than 
seen in the metal-rich sub-system of Galactic globular 
clusters, which has a width of about 0.3 dex. This metallic- 
ity width translates into a scatter of ~ 0.25 mag in V — K 



and ~ 0.15 mag in V — I. The red peak of the evolved 
intermediate-age clusters in NGC 4365 is at least twice as 
broad. However, both colour distributions, although sta- 
tistically insignificant, show hints for a substructure of 
the red peak which might consist of two or more red sub- 
populations (see the V — K colour histogram in Fig. ||). 

The intermediate-age sub-population of today's 
NGC 4365 globular clusters will mature to colours which 
are comparable to colours of the red (metal-rich) sub- 
population in M87 (Paper I). These broad red sub- 
populations are also found in other giant elliptical galax- 
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Table 4. Slopes oti, Pi as determined from the SSP models. The mean ay-K was calculated without the Vazdekis 
value. 



Model 



d/d(log t) /3» = <9/<9[Fe/H] 



parameter space 



Bruzual fe Chariot (2000) 



Kurth et al. (1999) 



Maraston (2001) 



Vazdekis (1999| ) 



V 

V ■ 

V ■ 
V- 

V - 
V- 

V ■ 

V - 



I 

K 

I 

K 

I 

K 

I 

K 



0.239 ± 0.029 
0.459 ±0.150 

0.218 ±0.043 
0.529 ±0.092 

0.260 ± 0.020 
0.554 ±0.017 

0.201 ± 0.052 
0.065 ±0.133 



0.236 ± 0.006 
0.964 ± 0.015 

0.259 ± 0.011 
0.981 ± 0.026 

0.207 ±0.011 
0.825 ± 0.020 

0.265 ± 0.006 
0.862 ± 0.030 



5 < t < 15 Gyr 
-0.7 <[Fe/H]< ±0.4 

5 < t < 15 Gyr 
-0.7 <[Fe/H]< ±0.4 

5 < t < 15 Gyr 
-0.35 <[Fe/H]< ±0.35 

5 < t < 15 Gyr 
-0.7 <[Fe/H]< ±0.4 



(mean) 



V -I 0.230 ± 0.026 
V-K 0.514 ±0.040 



0.242 ± 0.026 
0.908 ± 0.076 



5 < t < 15 Gyr 
-0.35 <[Fe/H]< ±0.35 



ies, such as NGC 4472 (e.g. peisler et al., 1996| ; [Puzia et 
al., 19^|), NGC 3585, NGC 4526, NGC 4649, NGC 5846 
(e-i 



Kundu 



Gebhardt & Kissler-Patig, 199£ ; Larsen et al. 
Whitmore, 2001a| ,|j|). 



2001 



3.5. Comparison with the Integrated Light of the Host 

Figure || shows the V — I and V — K colours of the in- 
tegrated light of NGC 3115 and NGC 4365 along with 
the colours of their globular cluster systems. The eff ective 

V — I colour was taken from |Buta fe Williams (1995| ). The 

V — K colour was measured by Frogel et al. (1978) in dif- 
ferent apertures which are not the effective apertures. We 
use the colour given by the aperture which is closest to 
the effective radius and calculate the expected change to 
an aperture of effective radius. All values are documented 
in Table |. 

From the comparison with SSP models both galaxies 



show diffuse stellar light colours which are consistent with 
a ~ 12 Gyr old Z > Zq metallicity stellar population. 
However, the observed colours are luminosity-weighted av- 
erages and can be reproduced by many different mixes of 
stellar populations with a wide range of ages and metal- 
licities. The major difference between the two galaxies ap- 
pears when the stellar light is compared with their glob- 
ular cluster systems (see Fig. ||). The globular clusters in 
NGC 3115 seem to have the same age as the galaxy body. 
The colours of the integrated light fall onto the metallicity 
sequence which is set up by the globular clusters. A few 
very red globular clusters share the colours of the stellar 
light. This clusters appear to be associated with the stellar 
disk of NGC 3115. 

In NGC 4365, the globular cluster system shows fea- 
tures that are not readily apparent in the diffuse stellar 
light. The galaxy's average stellar population is metal-rich 
and relatively old and the clusters span a wide range of 
metallicities and ages. 

The triaxial galaxy NGC 4365 ho sts a decoupled core 
which rotates around the minor-axis (Wagner et al., 1988 ; 
Bender fe Surma, 19921 ; purma fe Bender, 1995| ) while 



tailed abundance study found the innermost core stellar 
population being metal-rich 2.3Zq), relatively young 
7± 1.5 Gyr) and excessively enriched in a-elements 
(Surma & Bender, 1995). Davies et al. (2001) performed a 
subsequent integral-field spectroscopic study covering the 
galaxy out to ~ 3 kpc and found the same age (~ 14 Gyr) 
for the two distinct components, solar metallicity and a 
rather mild a-clcmcnt enhancement of the core. However, 
their Figure 2 shows that there is an age gradient in the 
very center (r < 1.6") where the innermost population is 
~ 12 Gyr old. The quite significant difference in age be- 
tween the former two studies can be attributed to the use 
of two different SSP models. 



The index measurements of both studies for H/3 and 
Mgl are consistent with each other. In particular, while 
purma fe Bender use Worthey (1992) models, Davies et 



al. (2001) use the models of (Vazdekis, 199£) which, in 



contrast to Worthey 's models, include late evolutionary 
phases, such as the HB (horizontal branch), AGB, EAGB, 
and PAGB. Since the HB can significant ly contribute to 
the H/3 flux ( |Maraston fe Thomas, 2000| ), which is com- 
monly used as an age-sensitive index, the predictions of 
a model which includes the HB will yield higher ages for 
same H/3 measurements as models without HB modeling. 
The difference in age between the previous two studies can 
therefore be fully understood from the choice of models. 



The latest Vazdekis99 models are also used in this 
study. Therefore, the spectroscopic and photometric pre- 
dictions for age and metallicity can directly be compared 
with each other. For NGC 4365, with the photometric 
VIK values of Table we obtain an age of 12 ± 4 Gyr 
and roughly solar metallicity. Considering that the pho- 
tometric values are an average over a large radial range, 
the predictions of photometric and spectroscopic parame- 
ters from the Vazdekis99 models are consistent with each 
other. 



the main body rolls around its major axis. The first de- 
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3.6. Comparison with the Milky Way and M31 
Globular Clusters 

We compare the colour-colour diagrams of NGC 3115 and 
NGC 4365 clusters with the ones of Milky Way and M31 
globular clusters. Figure || shows the comparison of the 
colour-colour diagrams. The globular cluster system of 
NGC 3115 compares well with the globular cluster system 
of either Local Group spiral. The globular cluster systems 
of NGC 3115, Milky Way, and M31 appear to be on aver- 
age old and roughly coeval. 

The few detected blue clusters in NGC 4365 show good 
agreement with the mean locus of metal-poor clusters in 
M31 and the Milky Way. The NGC 4365 globular clus- 
ter system extends also to much higher V — K at similar 
V — I colours than the Local Group systems, implying 
high metallicities and intermediate to young ages (see also 
Fig. |). 



4. Discussion 

4.1. NGC 4365, a Cooling- Down Galaxy 



of the intermediate-age globular cluster population to test 
how consistent these are with forming in the same event. 
To calculate the mass of the diffuse stellar component, we 



begin with the decoupled core. Davies et al. found that 



~ 12 Gyr tor the solar-mctallicity kinematically-decouplcd 
core (KDC) in NGC 4365 using the diagnostic p lot (H/3 
vs. [MgFe]) and SSP models of |Vazdekis (1999) ). From 
the comparison of our near-IR/optical colours with the 
same models, we find that ~ 40 — 80% of our sample are 
metal-rich [Z 0.5 Z®) globular clusters with ages be- 
tween ~ 2 — 8 Gyr. However, within our observational 
measurement errors the age distribution, derived from the 
linear interpolation of SSP models, is consistent with a 
single-burst of globular-cluster formation which took place 
~ 3 — 4 Gyr ago. Together with the metal-rich globular 
clusters a part of the stellar population in the KDC must 
have been formed from the same pre-enriched gas in a sig- 



nificant star-formation event. In fact, Davies et al. (2001) 
find that the ~ 2 Gyr younger luminosity-weighted age 
of the KDC (compared to the outer parts of the galaxy) 
can be accounted for by a contaminating, 5 Gyr old stel- 
lar population of solar metallicity. This young population 
would amount 6% of the mass inside r ^ 1.6" or ^ 160 
pc. 

It is tempting to associate the intermediate-age glob- 
ular cluster population with the kinematically decoupled 
core. Within the prediction accuracy of SSP models this 
young KDC population is consistent with being the dif- 
fuse counterpart of the stellar population in metal-rich 
intermediate-age globular clusters. Numerical simulations 
of gaseous mergers (see e.g. Barnes & Hcrnquist, 1992, for 
a review) predict the settling of the progenitors' gas into 
the inner parts of the merger remnant within ~ 10 8 years. 
If globular clusters and the central stellar population were 
in fact formed from the same gas, both stellar populations 
have to be almost coeval. 

We can compare the stellar mass contained in an 
intermediate-age component in NGC 4365 with the mass 



Davies! et al. (2001) find a luminosity-weighted age ot ^j" 



the intermediate-age component makes up 6% of the stel- 
lar mass of the decoupled core. Given the total stellar 
mass estimate for the core of M corc ^ (8.1 ± 2.5) • 10 9 Af© 
from Burma fc Bender (1995j ), the resulting mass in the 
core of the intermediate-age diffuse stellar component is 
(4.9 ±2.5) • 10 8 M Q . 

This mass is a lower limit to the total mass of the dif- 
fuse intermediate-age stellar population because it is un- 
likely that this component is found only in the very cen- 
tral region. In particular, the intermediate-age globular 
clusters are observed over a fairly extended region, and it 
would be difficult to arrange for the corresponding diffuse 
component not to have a similar extended spatial distri- 
bution. Therefore, in order to estimate the total stellar 
mass of the galaxy at intermediate ages, we need to esti- 
mate the fraction of the total galaxy mass coming from 
the intermediate age population. This fraction can not be 
greater than a few percent, or the extended intermediate- 
age population would have been detected by Davies et 
as they detected such a component in the central re- 
gion. Conversely, the mass fraction of the intermediate- 
age component must be at least a few-tenths of a percent 
based on comparing the intermediate-age stellar mass of 
(4.9 ± 2.5) • 10 8 Mq observed in the core to the total galaxy 
stellar mass of ~ 4.8 • lO n A/ ( [Bender et al., 1989Q . 

To bracket these possibilities, we carry out a represen- 
tative calculation with a 1% fraction of the overall stel- 
lar mass of the galaxy in an intermediate-age component, 
and show how the results are affected by varying this in 
either direction. This fraction gives a stellar mass over 
the whole galaxy in the intermediate-age component of 
-5.2-10 9 M Q . 

What is the corresponding mass currently in the 
intermediate-age GC system? After correcting our sample 
for incompleteness in luminosity and spatial coverage and 
with the assumption that ~ 40 — 80% of our sample globu- 
lar clusters have intermediate-ages, we find that the total 
number of intermediate-age globular clusters in NGC 4365 
is ~ 150 (within a factor of 2). To determine the mass in 
this system, we adopt the two-component mass function 
of McLaughlin & Pudritz (1996) over the mass range of 



10 — 10 Mq, appropriate for the GCs of giant elliptical 
galaxies. This gives a total mass of ~ 1.6 • 10 8 M Q in the 
globular cluster system. A similar result is found by simply 
counting up the estimated mass of the observed clusters 
and accounting for the spatial areas not covered by our 
observations. 

The efficiency of formation and survival to several Gyr 
of the intermediate-age globular cluster system is 3%. This 
follows directly from the derived mass of ~ 1.6- 10 8 M Q for 
the intermediate-age globular cluster system and the cor- 
responding ~ 5.2 • 10 9 M Q for the diffuse stellar component 
of the intermediate-age system. If the adopted mass frac- 
tion of the galaxy is too high, the formation and survival 
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1 2 3 4 12 3 4 

(V-K)„ fmagl (V-K), [magi 



Fig. 9. Comparison of globular clusters in NGC 3115 (left panel) and NGC 4365 (right panel) with globular clusters 
in the Milky Way and M31. Milky Way and M31 clusters are shown by open circles and open triangles. The stars show 
the colour of the diffuse stellar light (see Table |l|). As a guide to the eye, we plot a 12 Gyr iso-age line for different 
metallicities in the range —2.25 <[Fe/H]< +0.35 dex taken from Maraston (2001). The metallicity is increasing to 
redder colours and is indicated by solid squares. Note, that in the two Local Group spirals the high-metallicity range 
contains only few globular clusters which are, moreover, prone to have inaccurate reddening corrections. 



efficiency increases to rather high values of 10% or more. 
If the adopted mass fraction of the galaxy in intermediate- 
age component is a factor of a few too low, this number 
drops to about 1%. It is difficult to make the data consis- 
tent with a formation efficiency less than 1%. 

How do these formation and survival efficiencies com- 
pare to other observations? Observations of strongly star- 
bursting systems show that ~ 20% of the stars form in 



compact star clusters (e.g. Meurer et al 
et al 



1995|; |Whitmore 

19991; |Larsen fe Richtler, 1999|; |Zepf et al., 1999| ) 



This represents the peak formation efficiency of globular 
clusters. Studies of older mergers find th at a few percent 



of t he stars are in globular clusters (e.g. pchweizer et al. 
1996 |; [killer et al., 1997} |Zepf et al., 1999| ; |Goudfrooij et| 
al., 2001b). The best determined case for older globular 



cluster systems is probably the Galactic halo, in which 
about 1% of the stars are in globular clusters. A global 
aver age of all spheroida l galaxy systems is around 0.3% 
(e.g. McLaughlin, 1999 ). However, there are large varia- 
tions between galaxies, and within sub-populations in in- 
dividual galaxies. Dynamical evolution may significantly 
reduce the number of globular clusters observed at old ages 



compared to initial values at young ages (e.g. Gnedin & 
Ostrik^r, 19971; |Vcspcrini, 2000|) . 



The comparison of the formation and survival effi- 
ciency of the NGC 4365 intermediate-age globular clus- 
ter system to other systems suggests that NGC 4365 and 
older mergers have about 2% of their stars in globular 
clusters, with an uncertainty of a factor of a few. It is 



possible that the stellar mass attributed to the younger 
component should be restricted to that in the core alone, 
in which case the formation and survival efficiency of the 
NGC 4365 globular cluster system would be more like the 
peak efficiency seen in recent mergers with little subse- 
quent dynamical evolution. Underestimating the mass in 
the diffuse intermediate-age population would reduce the 
formation and survival efficiency by a factor of a few, but 
it is difficult to make the data consistent with globular 
cluster formation efficiencies less than a percent or so. 



We can also compare the metallicity of the 
intermediate-age globular cluster and diffuse stellar com- 
ponents. NGC 4365 shows super- metal- rich globular clus- 
ters (Z > 2 Zq) while the luminosity- weighted stellar light 
of the decoupled core reveals a roughly solar metallicity 
flSurma fe Bender, 1995[ [Davies et al., 200l| ). This is likely 
due to a mixing effect of the integrated light of stellar pop- 
ulations in the core. If stars and globular clusters were 
formed from the same pre-enriched material, super-solar 
metallicities are expected in the core, as well. The small 
fraction of light from super-metal-rich stars is likely to be 
washed out by insufficient resolution of the inner struc- 
tures. With a spatial r esolution of ~ 1.4" (rad ial size of the 
decoupled core ~ 8") Surma fc Bender (1995| ) report that 
the metallicity is at least 2.3Z in the center of NGC 4365. 
At face value, the core metallicity and the metallicity of 
the most metal- rich globular clusters are consistent. 
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4.2. NGC 3115, a Quiescent Galaxy 

It seems that the colour difference between the two ma- 
jor sub-populations of globular clusters in NGC 3115 is 
mainly due to a difference in metallicity, rather than in 
age. Within the errors this result is still consistent with 



the findings of Kundu & Whitmore (1998) who report, us- 



ing optical photometry only, that the metal-rich clusters 
are about ~ 3 Gyr younger than the metal-poor ones. 

In contrast to NGC 4365, NGC 3115 seems to have ex- 
perienced no major star-formation event in the last ~10 
Gyr. Both its globular cluster sub-populations are roughly 
that old or older. This resembles the situation in the glob- 
ular cluster system of the Milky Way, where the two sub- 
systems, corresponding to thick-disk/bulgc and halo popu- 
lations, differ mainly in metallicity ( Minniti, 1995| ; Cote et 
al., 20(j(], etc.) with a relatively small age difference where 
the metal-rich clusters may be slightly younger than their 



metal-poor counterparts (Rosenberg et al., 1999). 

Quantitatively, in NGC 3115 it remains to be seen 
whether or not the specific frequency Sn (number of glob- 
ular clusters per unit luminosity, see Harris & van den 



Bcrgh 1981) for the metal-poor stellar populations (i.e., 
the ratio of metal-poor globular clusters to metal-poor 
halo stars) is also roughly three times higher than for 
the metal-rich stellar populations, as it is observed in 
NGC 5128 and M31. 

Given the limited information about the halo stars, we 
conclude that the globular cluster system and the majority 
of halo stars in NGC 3115 most likely formed in two epochs 
and/or by two different mechanisms more than ~10 Gyr 
ago. Ever since then, the stellar populations in globular 
clusters and in the galaxy halo seem to have undergone a 
purely passive evolution. 

5. Summary 

Using accurate optical and near infrared photometry of 
the globular cluster systems in NGC 3115 and NGC 4365 
we found the following: 

— The large color baseline introduced by the new if-band 
data allowed us to find, for the first time, a very metal- 
rich intermediate-age population of globular clusters in 
NGC 4365. The clusters have ages between ~ 2 — 8 Gyr 
and metallicities ~ 0.5Zq—3Zq. They form a flattened 
population, which is aligned with the galaxy's major- 
axis, while the metal-poor clusters populate the halo. 
We emphasize that such intermediate-age populations 
can only be found from photometric data by obtaining 
deep X-band data in addition to high-quality optical 
imaging from HST/WFPC2. 

— After ~10 Gyr of passive evolution these globular clus- 
ters will show colours similar to super metal-rich glob- 
ular clusters in giant elliptical cluster galaxies, such as 
M87. 

— NGC 3115 hosts an old and bimodal globular clus- 
ter population. From comparison to four different SSP 



models we found the blue and the red globular cluster 
sub-populations being coeval to within the errors (±3 
Gyr). 

— We determine a metallicity difference between the red 
(metal-rich) and blue (metal-poor) globular cluster 
sub-population in NGC 3115 as A[Fe/H]= 1.0 ± 0.3 
dex. Thus, for this particular globular cluster system, 
metallicity is primarily responsible for the bimodality. 
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